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Abstract 

This  paper  addresses  the  modeling  of  strains  and  forces  generated  by  magnetostrictive  transducers  in 
response  to  applied  magnetic  fields.  The  magnetostrictive  effect  is  modeled  by  considering  both  the  ro¬ 
tation  of  magnetic  moments  in  response  to  the  field  and  the  elastic  vibrations  in  the  transducer.  The 
former  is  modeled  with  the  Jiles- Atherton  model  of  ferromagnetic  hysteresis  in  combination  with  a  quartic 
magnetostriction  law.  The  latter  is  modeled  through  force  balancing  which  yields  a  PDE  system  with 
magnetostrictive  inputs  and  boundary  conditions  given  by  the  specific  transducer  design.  The  solution  to 
this  system  provides  both  rod  displacements  and  forces.  The  calculated  forces  are  used  to  quantify  the 
magnetomechanical  effect  in  the  transducer  core,  i.e. ,  the  stress-induced  magnetization  changes.  This  is 
done  by  considering  a  ‘law  of  approach’  to  the  anhysteretic  magnetization.  The  resulting  model  provides 
a  representation  of  the  bidirectional  coupling  between  the  magnetic  and  elastic  states.  It  is  demonstrated 
that  the  model  accurately  characterizes  the  magnetic  hysteresis  in  the  material,  as  well  as  the  strains  and 
forces  output  by  the  transducer  under  conditions  typical  of  engineering  applications. 

Keywords:  Magnetostriction,  rare  earth-iron  alloys,  Terfenol-D,  magnetomechanical  model,  transducers 


1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

1999 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-1999  to  00-00-1999 

4.  TITLE  AND  SUBTITLE 

A  Coupled  Structural-Magnetic  Strain  and  Stress  Model  for 
Magnetostrictive  Transducers 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

North  Carolina  State  University, Center  for  Research  in  Scientific 
Computation, Raleigh, NC, 27695-8205 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

The  original  document  contains  color  images. 

14.  ABSTRACT 

see  report 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

ABSTRACT 

18.  NUMBER 

OF  PAGES 

24 

19a.  NAME  OF 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


1  Introduction 


We  address  the  modeling  of  the  magnetomechanical  behavior  of  transducer  systems  which  utilize  magne- 
tostrictive  materials  to  drive  structural  loads.  The  growing  interest  in  magnetostrictive  transducers  arises 
from  the  availability  of  highly  capable  magnetostrictive  materials,  such  as  the  rare  earth-iron  alloys,  which 
deliver  strains  in  the  10-3  range  and  forces  in  the  order  of  several  KN.  While  magnetostrictive  transducers 
provide  adequate  performance  at  the  low  signal  levels  where  their  behavior  is  quasilinear,  the  demand 
for  high  performance  transducers  often  dictates  that  they  be  driven  at  the  high  operating  regimes  where 
hysteresis  and  nonlinearities  are  intrinsic  to  magnetostrictive  behavior.  In  addition,  the  advantages  of 
magnetostrictive  materials  over  alternative  transducer  technologies  are  typically  realized  at  high  operat¬ 
ing  regimes.  This  motivates  the  development  of  models  that  accurately  characterize  the  hysteresis  and 
nonlinearities  intrinsic  to  magnetostrictive  transduction. 

Due  to  the  reciprocal  nature  of  magnetomechanical  transduction,  magnetostrictive  transducers  are  capa¬ 
ble  of  providing  not  only  actuation  but  sensing  capabilities  as  well.  One  crucial  aspect  of  the  bidirectional 
energy  exchange  between  the  magnetic  and  elastic  regimes  is  that  neither  transduction  mechanism  occurs 
independently,  but  rather  they  both  occur  simultaneously  during  transducer  operation.  In  the  direct  or 
magnetostrictive  effect,  the  action  of  a  magnetic  field  and  ensuing  magnetization  generates  strains  in  the 
magnetostrictive  material.  These  strains  are  in  turn  associated  with  a  stress  field  which  affects  the  magnetic 
state  by  virtue  of  the  inverse  or  magnetomechanical  effect.  This  issue  motivates  the  development  of  models 
capable  of  addressing  the  magnetomechanical  coupling  effects  intrinsic  to  magnetostrictive  transduction. 

The  close  connection  between  magnetostriction  and  the  magnetic  behavior  under  stress  has  long  been 
recognized,  and  extensive  experimental  evidence  on  the  magnetomechanical  effect  has  been  documented  [1- 
3] .  In  recent  years  there  has  been  renewed  interest  on  this  phenomenon  [4,  5]  because  of  its  relevance  for 
applications  such  as  non-destructive  evaluation  and  sensing.  In  the  specific  case  of  the  R-Fe2  (R  =  Tb,  Dy, 
Sm,  Ho)  alloys,  of  which  Terfenol-D  (Tbo.3Dyo.rFe1.9_2)  is  at  present  the  most  widely  known  commercially 
available  example,  it  has  been  demonstrated  that  desirable  properties  such  as  high  magnetostriction  and 
low  hysteresis  are  strongly  dependent  upon  the  operating  stress  [6-8]. 

Material  models  such  as  the  linear  piezomagnetic  equations  and  related  formulations  [9,  10]  are  widely 
used  to  characterize  magnetostrictive  transducer  performance.  The  most  common  formulation  for  these 
equations  is 

£  -  sH  a  +  d33  H  (1) 

B  =  d*33  a  +  ^H,  (2) 

in  which  e  is  the  strain,  sH  is  the  compliance  at  constant  field  H ,  d33  and  d33  are  the  magnetoelastic 
coupling  coefficients,  a  is  the  stress,  and  is  the  permeability  at  constant  stress.  It  is  emphasized  that 
this  model  is  in  essence  a  generalization  of  two  phenomenological  relationships,  namely  the  Hooke’s  law  for 
linearly  elastic  solids  e  —  sa  and  the  magnetic  constitutive  equation  B  =  j^iH.  The  total  magnetoelastic 
strain  e  given  by  equation  (1)  is  interpreted  as  the  superposition  of  the  elastic  or  passive  response  e  =  sa 
and  the  magnetostrictive  component  A  =  d33  H  associated  with  domain  processes  in  the  material.  In 
a  similar  fashion,  the  magnetic  induction  B  of  equation  (2)  is  interpreted  as  due  to  the  constant  stress 
magnetic  component  fia  H ,  and  a  term  due  to  magnetoelastic  interactions  d33a.  It  is  often  assumed  on 
the  basis  of  small  reversible  magnetostrictions  that  d33  =  d33,  which  suggests  that,  for  reversible  processes, 
a  large  magnetomechanical  effect  d33  =  (8B  /  da)jj  should  be  observed  in  materials  with  large  axial  strain 
coefficient  d33  =  ( de  /  dH)c r. 

The  linear  piezomagnetic  model  provides  adequate  characterization  of  magnetostrictive  performance 
at  the  low  operating  regimes  where  the  behavior  is  quasilinear  and  magnetic  hysteresis  is  negligible.  At 
high  operating  regimes,  however,  several  mechanisms  lead  to  hysteresis,  nonlinearities,  and  magnetoelastic 
coupling.  These  mechanisms  cannot  be  neglected  in  accurate  transducer  models. 
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This  paper  presents  a  nonlinear  and  hysteretic  magnetomechanical  model  for  the  strains  and  forces  gen¬ 
erated  by  magnetostrictive  transducers  in  response  to  applied  magnetic  fields.  The  formulation  presented 
here  extends  prior  modeling  capabilities  [11,  12],  which  include  magnetic  hysteresis  and  nonlinearities  as 
well  as  structural  vibrations,  by  including  the  sensing  or  magnetomechanical  effect  in  a  magnetostrictive 
rod  as  it  drives  a  transducer.  This  novel  modeling  approach  provides  an  improved  representation  of  the 
magnetomechanical  behavior  of  magnetostrictive  transducers  as  used  in  structural  applications,  and  it  pro¬ 
vides  the  framework  necessary  for  the  modeling  of  magnetostrictive  materials  when  utilized  as  sensors.  The 
model  is  illustrated  in  the  context  of  the  prototypical  Terfenol-D  transducer  depicted  in  Figure  1,  but  is 
not  limited  to  this  material  or  transducer  design.  This  design  provides  a  template  for  the  key  physical  com¬ 
ponents  required  to  fully  utilize  the  magnetostrictive  transducer  capabilities,  namely  the  magnetostrictive 
rod,  DC  magnetic  circuit,  AC  magnetic  circuit,  and  mechanical  prestress  mechanism. 

The  model  is  presented  in  three  stages.  In  the  first  stage,  discussed  in  Section  2,  we  consider  the  mag¬ 
netization  of  the  magnetostrictive  rod  under  an  externally  applied  magnetic  field  H  and  a  stress  field  a. 
The  field-induced  component  of  magnetization  is  quantified  with  the  mean  field  model  of  ferromagnetic 
hysteresis  originally  developed  by  Jiles  and  Atherton  [11-13].  The  stress-induced  component  of  magne¬ 
tization  is  modeled  with  a  law  of  approach  to  the  anhysteretic  magnetization  as  presented  in  [5].  The 
two  components  considered  together  provide  a  magnetization  model  based  on  the  energy  dissipated  when 
domain  walls  attach  to  and  unattach  from  inclusions  in  the  material.  Increasing  magnetic  fields  augment 
the  energy  dissipation  and  hence  the  amount  of  energy  lost  to  hysteresis.  Conversely,  increased  stresses 
help  to  overcome  domain  wall  pinning  and  hence  to  reduce  hysteresis. 

The  second  stage,  illustrated  in  Section  3,  involves  the  characterization  of  the  magnetostriction  A  pro¬ 
duced  when  the  magnetostrictive  rod  is  magnetized.  This  is  done  through  a  phenomenological  model 
consisting  of  an  even-terms  series  expansion.  While  A  includes  the  active  contribution  to  the  strain  arising 
from  the  rotation  of  magnetic  moments,  it  does  not  account  for  the  passive  or  material  response  of  the 
kind  found  in  ordinary  (i.e. ,  nonmagnetostrictive)  elastic  materials  and  modeled  by  sH  a  in  equation  (1). 

The  passive  effects  are  provided  in  Section  4  through  consideration  of  force  balancing  in  the  mag¬ 
netostrictive  rod,  in  the  form  of  a  PDE  equation  which  includes  the  intrinsic  magnetostriction,  system 
compliance,  internal  damping,  and  boundary  conditions  associated  with  the  mechanical  transducer  design. 
The  solution  to  this  PDE  provides  the  rod  displacements  and  corresponding  total  magnetoelastic  strain  e. 

The  complete  model  is  summarized  in  Section  5,  while  its  performance  is  evaluated  in  Section  6  by  means 
of  a  comparison  of  model  simulations  with  experimental  measurements.  Two  examples  are  provided.  In 
the  first,  the  accuracy  of  the  magnetization  and  strain  simulations  is  demonstrated  at  two  current  levels 
and  two  prestress  levels  for  a  transducer  designed  following  the  guidelines  of  Figure  1.  The  second  example 
demonstrates  the  ability  of  the  model  to  accurately  characterize  the  force  output  by  a  magnetostrictive 
rod  operated  under  mechanically  blocked  conditions  (e  =  0) . 


2  Magnetization  of  Magnetostrictive  Element 

Highly  magnetostrictive  materials  such  as  Terfenol-D  exhibit  substantial  deformations  when  magnetized. 
Changes  in  magnetization  are  due  to  the  application  of  magnetic  fields,  stresses,  or  thermal  energy.  At 
this  stage,  we  focus  on  changes  in  magnetization  and  subsequent  strains  produced  by  the  application  of 
magnetic  fields  and  stresses.  While  operating  temperature  has  a  strong  influence  on  the  performance  of 
magnetostrictive  materials  [16],  the  model  in  its  present  formulation  assumes  isothermal  behavior.  In  the 
presence  of  both  an  applied  magnetic  field  H  and  a  varying  stress  cr,  the  rate  of  change  of  magnetization 
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Figure  1:  Cross  section  of  the  prototypical  Terfenol-D  transducer  employed  for  model  development. 

with  time  dM  /  dt  is  dictated  by  the  expression 

dM  _  fdM\  dH  fdM\  da 
dt  V  9H )  dt  V  da  )  dt' 

which  motivates  treating  each  magnetization  term  independently.  In  the  first  term,  the  main  magnetoelastic 
component  is  the  differential  susceptibility  dM  /  dH ,  which  is  identified  from  the  ferromagnetic  hysteresis 
model.  The  time  rate  of  change  of  magnetic  held  dH  /  dt  is  readily  determined  since  the  magnetic  held 
input  to  the  transducer  is,  to  a  hrst  approximation,  known  in  advance.  In  the  second  term,  it  is  necessary 
to  quantify  both  the  magnetomechanical  effect  dM  /  da  and  the  time  rate  of  change  of  stress  da  /  dt,  since 
the  stress  depends  on  the  amount  of  strain  generated  by  the  transducer. 

2.1  Field-induced  magnetization:  differential  susceptibility 

The  magnetization  of  a  ferromagnetic  material  in  response  to  applied  magnetic  helds  can  be  explained  by 
considering  two  related  mechanisms  [15-18].  The  hrst  mechanism  is  that  domain  walls  (the  transition  layers 
between  highly  aligned  regions  termed  magnetic  domains)  move  under  the  influence  of  the  magnetic  held 
in  such  a  way  that  favorably  oriented  domains  grow  at  the  expense  of  unfavorably  oriented  domains.  The 
wall  thickness  is  determined  by  a  balance  between  the  anisotropy,  which  tends  to  make  the  walls  thinner, 
and  the  Weiss-type  interaction  coupling  between  atomic  magnetic  moments,  which  tends  to  make  the  walls 
thicker.  The  second  mechanism  involves  the  rotation  of  magnetic  moments  within  domains  towards  the 
held  direction. 

In  the  theory  originally  proposed  by  .files  and  Atherton  [11-13,19]  and  implemented  here,  the  magnetiza¬ 
tion  is  quantified  from  the  difference  between  the  total  magnetic  energy  available  to  magnetize  the  material 
and  the  energy  lost  to  domain  wall  pinning.  The  former  energy  state  is  achieved  in  the  anhysteretic  condi¬ 
tion,  which  has  a  magnetization  value  Man.  The  energy  lost  to  pinning  is  modeled  assuming  a  friction-type 
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mechanism  which  opposes  changes  in  magnetization.  To  quantify  these  energies,  minimization  of  a  suitable 
thermodynamic  potential  A  is  used  to  identify  first  the  effective  field  He  acting  on  the  material, 

He  —  H  +  a  M  +  Ha)  (4) 


where  H  is  the  applied  field,  a  M  is  the  Weiss  interaction  field  responsible  for  the  alignment  of  neighboring 


■  a  e 


jdM  is  the  field  due  to  magnetoelastic 


magnetic  moments  within  domains,  and  Ha  =  l/f-io  Q 
interactions. 

The  effective  field  calculated  from  expression  (4)  is  then  used  to  compute  the  anhysteretic  magnetization, 
which  is  quantified  using  the  Langevin  function  C(z )  =  coth(z)  —  1/z,  with  —  1  <  C(z)  <  1, 


Man  =  MsC{He/a)} 


(•5) 


in  which  Ms  is  the  saturation  magnetization  and  the  constant  a,  representing  the  effective  domain  density, 
is  treated  as  a  parameter  to  be  estimated  through  a  least  squares  fit  to  the  data.  It  is  emphasized  that 
while  the  Langevin  function  provides  satisfactory  fits  of  the  ferromagnetic  anhysteretic  magnetization,  its 
applicability  to  ferromagnets  should  be  considered  as  semi-phenomenological  in  nature.  The  Langevin 
function  was  originally  developed  for  the  paramagnetic  state,  in  which  all  possible  orientations  have  the 
same  probability  and  the  material  can  be  considered  isotropic.  In  addition  to  being  intrinsically  anisotropic, 
the  motion  of  domain  walls  in  ferromagnetic  materials  is  impeded  by  imperfections  or  pinning  sites.  These 
pinning  sites,  which  are  attributable  to  crystallographic  defects  and  to  the  presence  of  dendritic  twin 
boundaries  in  the  case  of  Terfenol-D,  form  energy  wells  that  are  energetically  favorable  for  domain  wall 
attachment. 

The  effect  of  pinning  sites  on  domain  wall  motion  under  constant  stress  conditions  has  been  assessed 
through  consideration  of  reversible  Mrev  and  irreversible  M{rr  components  of  the  magnetization.  For  low 
magnetic  field  intensities  about  an  equilibrium  level,  the  domain  walls  bend  reversibly  while  remaining 
attached  to  pinning  sites.  As  the  applied  magnetic  field  is  increased,  the  domain  walls  achieve  sufficient 
energy  to  break  free  from  pinning  sites  while  moving  up  the  energy  well  where  they  were  originally  located, 
and  attach  to  remote  sites  where  the  energy  configuration  is  favorable.  The  energy  lost  to  domain  wall 
pinning  manifests  itself  as  hysteresis  in  the  magnetization.  This  means  that  under  the  assumption  of  no 
other  loss  mechanisms,  the  hysteresis  loss  per  unit  volume  and  per  cycle  can  be  quantified  experimentally 
from  the  area  enclosed  by  the  M-H  loop. 

Energy  balancing  is  used  to  derive  a  differential  equation  for  irreversible  changes  in  magnetization,  which 
can  be  shown  to  be  [15] 


Mirr  =  Ma 


k  5 


d.\ijrr 

dHe 


(6) 


where  the  constant  k  quantifies  the  energy  needed  to  break  pinning  sites,  and  S  has  the  value  +1  when 

dH  /  dt  >  0  and  -1  when  dH  /  dt  <  0  to  ensure  that  the  energy  lost  to  pinning  always  opposes  magnetization 

changes.  Applying  the  chain  rule,  expression  (6)  can  be  modified  to  give  the  differential  irreversible 
susceptibility, 

dMirr  _  Man  -  Mirr  9He 

dH  ~  Sk  dH'  U 


Recognizing  that  in  this  case  the  effective  field  given  by  equation  (4)  should  be  defined  in  terms  of  the 
irreversible  magnetization, 


H 


H  +  a  Mirr  + 


3  d  (a  s) 

2  fio  dMirr  ’ 
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the  partial  derivative  of  the  effective  field  with  respect  to  field  takes  the  form 


dHe 

dH 


1  -(“  O  (  .1  I  ,  O  ) 


dMirr 

oh 


where  a(Mirr,o)  is  an  unitless  effective  coupling  term  which  is  defined  as 


a{Mirr ,  a)  =  a  + 


3  d2(crs) 

2^  dMfr  ' 


(8) 


It  is  noted  that  the  final  form  of  a  depends  upon  the  characterization  of  £  and  a. 

Direct  substitution  of  (8)  into  (7)  yields  the  differential  equation  for  the  irreversible  susceptibility, 


dMirr  =  _ Man  ~  Mirr _ 

OH  5k  -  a (Mirr,  ct) (Man  -  Mtrr) ' 


(9) 


It  is  noted  that  equation  (9)  can  yield  nonphysical  solutions  when  dH  is  reversed  near  saturation.  Specifi¬ 
cally,  when  the  magnetization  M8rr  is  below  the  anhysteretic  Man  for  increasing  field,  or  when  the  magne¬ 
tization  Mirr  is  above  the  anhysteretic  Man  for  decreasing  field,  direct  solution  of  (7)  leads  to  a  negative, 
and  hence  nonphysical,  differential  susceptibility.  A  mathematical  strategy  has  been  devised  in  [22]  that 
produces  a  more  physically  consistent  expression 

dMjrr  _  c  _ Man  ~  Mirr _ 

dH  8  5  k  -  a (Mirr,  ct)  (Man  -  Mirr) 

where 

1,  {dH  /  dt  >  0  and  M  >  Man}  or  {dH  /  dt  <  0  andM  <  Man} 

0,  otherwise. 

Numerical  integration  of  differential  equation  (10)  gives  the  irreversible  magnetization  Mjrr. 

It  has  been  hypothesized  in  [15]  that  the  reversible  component  of  magnetization  reduces  the  difference 
between  the  prevailing  irreversible  magnetization  .V/,,T  and  the  anhysteretic  magnetization  Man  at  the 
same  field.  This  can  be  modeled  mathematically  with  the  expression 


(10) 


Mrev  =  c(Man- Mlrr),  (11) 

where  the  coefficient  of  proportionality  c  quantifies  the  amount  of  reversible  domain  wall  bulging.  The 
value  of  c  is  determined  experimentally  from  the  ratio  of  the  initial  and  anhysteretic  susceptibilities  [22] 
or  through  a  fit  to  experimental  data.  Differentiation  of  equation  (11)  gives  the  reversible  differential 
susceptibility 

dMrev  _  f  dMan  dMirr  \ 

dH  ~c  v  dH  dH  J 

Finally,  the  irreversible  and  reversible  terms  given  by  expressions  (10)  and  (12),  are  added  to  give  the 
total  susceptibility 

OM  _  _ Man  ~  Mjrr _  dMan 

dH  1  C)  4  5  k  -  a (Mirr,  <r)  (Man  -  Mtrr)  +  °  dH  ’  j 

which  after  numerical  integration  gives  the  total  magnetization  arising  from  the  application  of  a  magnetic 
field.  Alternatively,  M  can  be  computed  by  solving  equation  (10)  and  then  adding  Mrev  from  equation  (11). 


(12) 
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2.2  Stress-induced  magnetization  changes:  magnetomechanical  effect 

We  now  consider  the  contribution  of  stress  to  the  total  magnetization,  or  magnetomechanical  effect 
DM  /da.  A  unifying  description  of  the  changes  in  magnetization  due  to  the  action  of  stress  has  been 
recently  developed  [4,  5,  13].  In  the  theory  presented  by  Jiles  [5,  14]  and  implemented  here,  the  main 
mechanism  governing  the  magnetomechanical  effect  is  the  unpinning  of  domain  walls  produced  when  a 
stress  is  applied.  On  the  basis  of  the  key  model  assumption  that  hysteresis  is  originated  mainly  from 
domain  wall  pinning,  the  freeing  of  domain  walls  from  their  pinning  sites  must  cause  the  magnetization  to 
change  in  such  a  way  as  to  approach  the  anhysteretic. 

Experimental  measurements  reported  in  [3,  4,  14,  23]  clearly  indicate  that  both  the  magnitude  and 
direction  of  stress-induced  magnetization  changes  are  profoundly  influenced  by  the  magnetic  history  of  the 
specimen.  To  illustrate,  Figure  2  shows  a  schematic  representation  of  the  approach  to  the  anhysteretic. 
The  sample  was  first  magnetized  to  saturation,  and  the  field  was  then  removed  so  that  the  magnetization 
lies  close  to  positive  or  negative  remanence,  as  respectively  denoted  by  points  A  and  B  in  Figure  2(a). 
The  magnetization  changes  AM  were  then  computed  for  varying  compressive  stresses  while  the  field  was 
held  constant.  The  magnetization  changes  are  illustrated  in  Figure  2(b).  ft  is  noted  that  the  magnitude 
and  direction  of  AM  depend  on  whether  the  initial  magnetization  lies  above  or  below  the  anhysteretic. 
Furthermore,  the  magnetization  at  point  X  changes  slightly  upon  removal  of  the  stress. 

The  change  in  magnetization  exhibited  by  point  X  upon  removal  of  the  stress  is  explained  by  a  second, 
related  effect,  that  is  the  anhysteretic  curve  itself  varies  under  the  action  of  stress.  Application  of  tensile 
stress  to  a  material  with  positive  magnetostriction  coefficient  (such  as  Terfenol-D)  produces  an  increase  in 
both  the  slope  of  the  M-H  loop  and  the  remanent  magnetization  value,  while  compressive  stress  produces 
a  shearing  of  the  M-H  loop.  This  effect  is  illustrated  in  Figure  3,  where  the  influence  of  stress  on  both  the 
anhysteretic  and  total  magnetization  is  shown  for  the  stress  levels  +10,  0,  -10  and  -20  MPa. 

The  net  result  of  both  effects  considered  in  combination  is  that  a  monotonically  increasing  stress  causes 
the  magnetization  value  of  X  to  be  approached  in  the  fashion  depicted  in  Figure  2,  while  the  magnetization 
value  of  X  itself  varies  as  shown  in  Figure  3.  The  magnetomechanical  effect  model  utilized  here  is  then 
formulated  in  the  context  of:  (i)  the  effect  of  stress  on  the  anhysteretic  magnetization  and  (ii)  a  law  of 
approach  to  the  anhysteretic  magnetization  upon  application  of  stress. 

2.2.1  Anhysteretic  Magnetization 

The  anhysteretic  magnetization  Man  given  by  expression  (5)  should  be  interpreted  as  a  local  function 
of  both  H  and  the  particular  value  of  the  prevailing  magnetization  M,  with  both  operating  through 
the  effective  field  He  given  by  (4).  For  instance,  at  point  A  in  Figure  4,  the  material  is  in  a  state  of 
equilibrium  in  the  presence  of  zero  applied  field,  due  to  the  coupling  between  magnetic  moments  a  M 
and  the  magnetoelastic  interactions  H a  in  the  material.  Application  of  a  small  magnetic  field  perturbs 
this  equilibrium  slightly,  giving  a  new  value  of  Man.  Applying  small  field  increments  successively,  the 
two-valued  local  anhysteretic  curve  shown  in  Figure  4  is  obtained.  Note,  however,  that  global  equilibrium 
is  never  achieved  in  this  case. 

A  closer  look  at  the  Langevin  function  (5)  reveals  that  it  is  possible  to  find  a  solution  Man  which,  for 
a  given  value  of  H  and  <r,  satisfies  the  equation  identically, 

Man  =  MsC[He(Man)/a\.  (14) 

The  solution  Man  to  the  Langevin  function  is  unique,  and  may  be  interpreted  as  a  state  of  equilibrium 
associated  with  the  minimization  of  the  total  energy.  In  consequence,  the  locus  of  points  obtained  by 
mapping  all  possible  values  of  H  for  a  given  stress  state  is  a  single-valued  function,  as  illustrated  in 
Figure  4.  We  term  this  function  the  ‘global’  anhysteretic  magnetization.  The  global  anhysteretic  Man 
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(a)  (b) 


Figure  2:  Representation  of  the  approach  to  the  anhysteretic  under  stress,  after  [4];  -  total 

magnetization,  —  •  —  anhysteretic  magnetization,  (a)  Arrows  indicate  approach  in  the  M-H  plane 
from  positions  A  and  B  above  and  below  the  anhysteretic,  for  fixed  H .  In  either  case,  M  moves 
towards  point  X  on  the  anhysteretic  curve,  (b)  Trajectory  of  magnetization  change  A M/Ms  upon 
application  and  further  removal  of  a  compressive  stress,  starting  at  point  A  (A  — >  X),  and  at  point 
B  (B  X). 
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Figure  3:  Model  simulations  representing  the  effect  of  stress  on  the  total  magnetization  ( - ), 

and  on  the  global  anhysteretic  magnetization  (—  •  — ),  for  +10,  0,  -10,  and  -20  MPa.  A  positive 
magnetostriction  coefficient  is  assumed  [14], 


is  the  anhysteretic  used  to  describe  the  law  of  approach,  on  the  basis  of  the  observation  [4]  that  the 
application  of  stress  causes  the  prevailing  magnetization  to  approach  the  single-valued  anhysteretic  state 
obtained  in  experiments  by  superimposing  a  decaying  AC  field  on  top  of  a  fixed  DC  field.  It  is  noted  that 
in  the  absence  of  magnetic  coupling  (o  =  0)  and  constant  magnetoelastic  interactions  ( Ha  ^  Ha ( M ) ) .  the 
local  and  global  anhysteretic  curves  coincide. 

The  global  anhysteretic  depends  significantly  on  the  stress  state  in  the  material.  It  has  been  even 
suggested  [24]  that  the  global  anhysteretic  defines  not  a  curve,  but  a  surface  dependent  on  H  and  a.  This 
dual  dependency  is  described  in  terms  of  the  effective  field  H, ,  by  recognizing  [5]  that  the  global  anhysteretic 
magnetization  under  a  field  H  and  stress  a  is  equivalent  to  the  global  anhysteretic  magnetization  under 
field  He  and  zero  stress,  M an(H,  a)  —  .V/ ,  0). 

In  order  to  compute  the  global  anhysteretic,  the  effective  field  given  in  equation  (4)  is  written  in  the 
form 


He  —  H  +  a(M,  a)  M, 

where  the  effective  coupling  parameter  a(M,  a)  is 

,ir  ,  3  1  d(ae) 

v  '  2  no  M  dM 

The  global  anhysteretic  is  then  given  by  the  following  expression, 

'  H  +  a (Mam  a)  M an  \  ( 


Man  =  M, 


coth 


a  J  \H  +  a(Man,  a)  Man  ^ 

which  may  be  solved  numerically  using  an  iterative  technique  such  as  the  Newton-Raphson  method. 


Figure  4:  Model  simulations  showing  the  relative  magnetization  M/MS)  relative  local  anhysteretic 
Man(H,  M,  a) /Ms  =  cot h [( H  +  5(M,  a)M)/a]  —  \a/(H  +  5(M,  a)M)],  and  relative  global  anhysteretic 
Man{H,a)/Ms  =  coth[(i?  +  a(Man,  a) Man)/ a]  ~  W/{H  +  a(Man,a)Man)],  as  a  function  of  H. 
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2.2.2  Law  of  approach 


It  has  been  observed  experimentally  that  the  direction  in  which  the  magnetization  changes  with  applied 
stress  is  independent  of  the  sign  of  the  stress,  for  small  stresses  and  when  the  magnetization  is  sufficiently 
distant  from  the  anhysteretic.  It  is  then  inferred  that  the  direction  of  change  is  dependent  not  on  the  stress 
itself,  but  on  a  quantity  which  is  independent  of  the  sign  of  the  stress.  In  this  context,  Jiles  hypothesized  [5] 
that  this  quantity  is  the  elastic  energy  per  unit  volume,  W  =  a2/(2E))  which  is  clearly  independent  of 
the  sign  of  a.  The  ‘law  of  approach’  to  the  anhysteretic  condition  is  then  formulated  as  follows:  the 
rate  of  change  of  magnetization  with  elastic  energy  is  proportional  to  the  displacement  of  the  prevailing 
magnetization  from  the  global  anhysteretic  magnetization.  The  concept  of  the  law  of  approach  is  now 
applied  to  the  stress-induced  magnetization  of  a  magnetostrictive  material. 

As  before,  this  may  be  posed  via  irreversible  and  reversible  components  of  the  magnetization.  It  is  noted 
that  to  a  first  approximation,  the  application  of  stress  produces  irreversible  magnetization  changes  since 
AM  arising  from  stress  unloading  is  negligible.  Thus,  it  is  reasonable  to  formulate  the  law  of  approach  in 
terms  of  the  irreversible  magnetization  Mvrr . 

=  I  (M„„  -  M,rr),  (15) 

where  4  is  a  coefficient  with  dimensions  of  energy  per  unit  volume  that  needs  to  be  identified  for  magne¬ 
tostrictive  materials.  Making  the  substitution  dW  /  da  =  a  j  E  in  (15),  along  with  application  of  the  chain 
rule,  permits  writing  the  irreversible  component  of  the  magnetomechanical  effect  as  follows, 

dM;rr  a  , — 

(Man  ~  Mirr).  (16) 

A  similar  argument  to  that  used  in  the  field-induced  case  yields  an  expression  for  the  reversible  compo¬ 
nent, 


dMrev  _  ( dMan  0Mrrr  \ 

da  ~  C  ^  da  ~  da  )' 


(17) 


with  c  the  same  coefficient  as  that  defined  in  equation  (11)  because  the  energy  available  for  domain  wall 
bulging  should  be  independent  of  the  mechanism  that  produces  the  bulging,  which  can  be  either  field-  or 
stress-induced. 

Summing  the  irreversible  and  reversible  contributions  leads  to 


dM 

Ih 


(1  —  c)  T77  (Man  ~  Mjrr)  +  C  ■ 

E  4 


dMn 


da 


(18) 


which  after  numerical  integration  gives  the  magnetization  arising  from  the  application  of  stress.  As  be¬ 
fore,  M  may  be  alternatively  computed  by  solving  M{rr  from  equation  (16)  and  then  adding  Mrev  from 
equation  (17). 


3  Active  Component  of  Strain 

It  is  ultimately  necessary  to  quantify  the  strains  generated  by  the  magnetostrictive  material  when  a  mag¬ 
netic  field  is  applied.  For  this  purpose,  we  consider  first  the  deformations  occurring  in  the  crystal  lattice 
when  the  domain  configuration  changes.  Several  models  exist  for  quantifying  these  deformations,  includ¬ 
ing  the  quadratic  law  for  domain  magnetization  rotation  discussed  in  [19,  20],  energy  or  thermodynamic 
formulations  [23-25],  elastomagnetic  models  [26-30],  micro  magnetic  theories  [32],  magnetization  rotation 
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analysis  [7],  and  empirical  models  [5].  At  low  to  moderate  operating  levels,  or  when  material  stresses  are 
invariant,  these  deformations  dominate  over  other  material  elastic  dynamics.  Under  such  circumstances,  it 
is  theoretically  possible  to  quantify  the  bulk  magnetostriction  upon  knowledge  of  the  domain  configuration 
and  the  magnetostriction  along  easy  crystallographic  axes.  In  the  case  of  Terfenol-D,  nominal  values  for 
the  latter  are  Am  =  1600  X  10-6  and  Aioo  =  90  X  10-6,  and  As  m  1000  X  10-6.  In  practical  terms,  however, 
the  domain  configuration  cannot  be  known  apriori. 

To  motivate  the  approach  followed  in  this  work,  we  consider  the  particular  case  when  the  magnetic  field 
is  applied  perpendicular  to  the  axis  in  which  the  magnetic  moments  have  been  aligned  by  application  of 
sufficiently  large  compression  in  the  case  of  a  poly  crystalline  material  such  as  Terfenol-D,  or  perpendicular 
to  the  easy  crystallographic  axis  in  a  single  crystal  with  uniaxial  anisotropy.  In  either  case  domain  rotation 
is  the  prevailing  magnetization  mechanism,  and  the  magnetostriction  along  the  field  direction  is  given 
by  [19] 


m)  =  \  a. 


(19) 


which  predicts  a  quadratic  relation  between  A  and  M .  Expression  (19)  is  a  single- valued  functional,  while 
extensive  experimental  evidence  demonstrates  that  the  A -M  relationship  exhibits  some  degree  of  hysteresis. 
For  transducer  modeling  purposes,  it  is  feasible  to  utilize  a  single  valued  A -M  functional  to  model  the  overall 
shape,  and  to  let  M  provide  the  hysteresis  through  the  hysteretic  mechanisms  in  M-H .  This  approach  has 
proven  effective  in  previous  investigations  [33]. 

Even  though  (19)  is  consistent  with  the  physical  phenomena  occurring  in  the  above  mentioned  cases,  it 
is  not  sufficiently  general  in  cases  when  domain  wall  motion  is  significant  such  as  when  the  operating  stress 
acting  on  the  Terfenol-D  material  is  not  extreme.  In  order  to  provide  a  more  general  magnetostriction 
model,  we  consider  a  series  expansion  symmetric  about  M  =  0, 


A  (M)  =  ^7i'M24, 
*= o 


in  which  the  coefficients  7 j  need  to  be  identified  from  experimental  data.  It  is  noted  that  quadratic 
relation  (19)  is  achieved  for  i  =  1  with  70  =  0  and  71  =  (3  AS)/(2MS2).  For  implementation  purposes,  we 
consider  in  this  study  a  quartic  law  in  which  the  series  is  truncated  after  i  =  2, 

A(M)  =  71  M2  +  72  M4.  (20) 


Note: 

A  convenient  method  for  identifying  71  and  72  is  by  solving  simultaneously  the  constraints  at  magnetic 
saturation  Ms  and  at  the  inflection  point  Mo, 

f  A(Ms)  =  7i  Mg  +  72  Mg  = 

1  d2X 

{  qm2(M 0)  =  2 7l  +  12 72  Mq  =  0, 

in  which  (3/2)  As  is  the  value  of  the  magnetostriction  at  saturation  magnetization. 

4  Passive  Component  and  Total  Strain 

The  magnetostriction  A  given  by  (20)  quantifies  the  reorientation  of  magnetic  moments  towards  the  di¬ 
rection  of  applied  field.  It  was  shown  in  [12]  that  this  magnetostriction  is  analogous  to  the  term  d^H  in 
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linear  models,  but  is  inherently  nonlinear  and  hysteretic  through  both  the  magnetization  M  and  the  quar- 
tic  relation  between  M  and  A.  It  ignores,  however,  the  elastic  properties  of  the  magnetostrictive  material 
as  it  vibrates,  as  represented  in  the  linear  models  by  sHa.  In  this  section,  a  PDE  system  is  formulated 
which  models  the  elastic  response  of  the  magnetostrictive  material  and  relevant  transducer  components 
located  in  the  load  path.  The  input  to  this  PDE  is  formulated  through  the  magnetostriction  A,  which 
constitutes  an  ‘internal  force’  driving  the  vibrations  of  the  transducer.  The  solution  to  the  PDE  is  the 
longitudinal  displacements  u(t ,  x )  relative  to  the  prestressed  position.  Additional  details  regarding  this 
PDE  formulation  are  provided  in  [12]. 

The  structural  dynamics  are  modeled  through  consideration  of  the  magnetostrictive  rod  (5),  prestress 
bolt  (6),  spring  (7),  and  mass  load  (8)  for  the  transducer  in  Figure  1.  The  prestress  bolt  provides  a  stress 
(To  <  0  by  compressing  the  magnetostrictive  rod  against  the  spring,  modeled  by  a  linear  spring  k /  and 
dashpot  cl.  The  rod  is  assumed  to  have  length  L ,  cross  sectional  area  A,  and  longitudinal  coordinate  x. 
The  material  density  is  p.  the  elastic  modulus  is  E ,  the  internal  damping  is  cp ,  and  the  external  load  is 
modeled  by  a  point  mass  rn /  , .  It  should  be  noted  that  parameter  E  lies  between  the  elastic  modulus  at 
constant  H ,  EH ,  and  at  constant  B ,  EB .  Since  EH  and  EB  depend  upon  the  field  intensity  [34],  so  does 
E.  However,  for  simplicity  E  is  treated  as  a  nominal  or  operational  material  stiffness. 

As  detailed  in  [11,  12],  the  total  stress  at  cross  sections  x  in  the  rod  is  given  by 

du  d2  u 

a(t.,x)  =  E—  (t,  x)  +  cD  ~  EX(t,x)  +  a0,  (21) 

where  the  terms  on  the  right  hand  side  represent  respectively  the  linear  elasticity  at  small  displacements, 
Kelvin- Voigt  damping,  magnetostriction-derived  stress,  and  prestress.  Force  balancing  then  yields  the 
dynamic  model  for  the  longitudinal  displacements  and  strains. 

For  implementation  purposes,  the  model  is  formulated  in  weak  or  variational  form  by  multiplying  the 
strong  form  by  test  functions  0  followed  by  integration  throughout  the  length  of  the  rod.  This  reduces 
the  smoothness  requirements  on  the  finite  element  basis  since  displacements  and  test  functions  need  be 
differentiated  only  once  compared  to  the  second  derivatives  present  in  the  strong  form.  The  space  of  test 
functions  is  V  =  ifjl  ( 0,  L)  =  {(f>  €  H1( 0,  L )  |  <^(0)  =  0},  so  that  for  all  <f>{x)  G  V, 

fL  .  d2u  .  fL  T  du  d2u  ,  ,  w  1  dcj) , 

j0  PA-g^{t^X)HX)dx  =  -  Jo  EAg^(^x)  +  cDAg^{t,x)-EAX(t,x)  —  (*)  rf* 

Dv, 

-  kLii(t,L)  +  cL—(t,L)  +  m,L—(t,L)  (22) 

The  solution  u(t,x)  to  expression  (22)  defines  the  longitudinal  displacements  about  the  prestressed 
position.  Once  the  displacements  are  computed,  the  strains  are  evaluated  by  taking  derivatives  with 
respect  to  position,  e(t,  x )  =  du  /  dx (t,  x),  and  the  material  stresses  a(t,  x)  are  calculated  directly  from  (21). 
Note  that  the  stress  at  the  rod  end  a(t,  /.)  may  be  equivalently  calculated  from  the  boundary  condition 
ait^  L )  =  ji~kL  «(*,  L)  ~  cl  | f  (T  L)  -  mL  | ^{t,  L)]. 

5  Transducer  Model  Summary 

In  the  presence  of  both  an  applied  magnetic  field  H  and  stress  a,  the  total  magnetization  is  dictated 
by  the  superposition  of  the  field-  and  stress-dependent  components  given  by  equations  (13)  and  (18). 
For  implementation  purposes,  however,  we  consider  the  alternative  approach  in  which  the  irreversible 
magnetization  Mirr(t,x )  is  computed  first  and  added  to  the  reversible  component  ,V/,.f ,  (/.  ./•)  to  give  the 
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total  magnetization  M(t,x).  The  superposition  of  expressions  (10)  and  (16)  leads  to  the  differential 
equation  for  the  time  rate  of  change  of  irreversible  magnetization: 

dMjrr  _  f  _ Man(t.,x )  -  Mirr(t,x) _ \  dH 

dt  ’  '  1  8  k  —  a (Mirr,a)  [Man(t,  x)  -  Mirr(t,x)]  J  dt  ’ 

+  |  —JTJ—  \Man(t,  x)  -  Mirr(t,x)  j^(t,  x).  (23) 

To  characterize  dH  /  dt,  it  is  necessary  to  quantify  first  the  quasistatic  field  H(t,x )  generated  by 
the  solenoid  when  a  current  /(f)  circulates  through  it.  To  this  end,  it  is  often  assumed  that  //(f)  = 
(No.  turns  /length)  /(f).  However,  this  model  is  only  valid  in  the  idealized  situation  of  a  lossless,  infinitely 
long  solenoid  in  a  lossless  magnetic  circuit.  Experimental  evidence  on  research  transducers  indicates  that 
H  =  n  I  is  highly  inaccurate  because  it  neglects  solenoid  end  effects,  demagnetizing  factors,  ohmic  losses, 
and  flux  leakage.  One  possible  approach  consists  of  identifying  H-I  by  solving  numerically  Ampere’s  law 
or  the  Biot  Savart  law,  using  for  instance  finite  element  methods.  For  purposes  of  implementing  the  cou¬ 
pled  magnetomechanical  model,  an  experimental  approach  was  followed  which  consisted  of  experimentally 
determining  the  magnetic  circuit  behavior  via  a  position-dependent  filter  $(a:)  which  accounts  for  the 
above-mentioned  effects.  Thus,  the  time  and  spatial  dependencies  of  the  field  are  formulated  with  the 
expression 

H(t,  x)  =  I(t) 'lt(x).  (24) 

Upon  substitution  of  equation  (24)  into  (23),  the  final  form  for  the  rate  of  change  of  irreversible  mag¬ 
netization  is  obtained, 

dMirr  ,  .  _  j  Man  (f,  x)  Mj>r(f,  x)  1  dl,  \T/ f  £*  i 

dt  ’  \  8  k  —  a(Mirr,  a)  [Man{t,  x)  —  Mirr(t,  x)]j  dt 

+  '  \Man(t,x)  -  Mtrr(t,x)  }^-(f,a0,  (25) 

which  relates  the  irreversible  magnetization  with  the  current  applied  to  the  transducer  and  the  varying 
stress  arising  during  operation. 

The  reversible  magnetization  is  calculated  directly  upon  integration  and  subsequent  superposition  of 
equations  (11)  and  (17), 

Mrev(t,  x)  =  c  [Man(f,  x)  +  Man(t,  x)  -  cMirr(t,  x).  (26) 

The  total  magnetization  under  the  application  of  both  a  field  and  a  stress  is  then  given  by 

M(f,  x)  -  (1  —  c )  Mirr(t,  x)  +  c  \[Man{t1  x)  +  Ma„(f,  x )  , 

which  includes  the  irreversible  and  reversible  components  given  by  equations  (25)  and  (26),  respectively, 
ft  should  be  noted  that  in  the  case  of  constant  stress  (da  /  dt  =  0)  or  constant  field  ( dl  /  dt  =  0) ,  the 
expression  reduces  to  the  individual  components  characterized  by  expressions  (13)  and  (18). 

After  the  magnetization  M(t,x )  arising  from  the  application  of  H(t.,x )  and  a(t.,x )  has  been  identified, 
the  active  component  of  strain  is  computed  from  equation  (20), 

A [M(t,  x)\  =  7!  M2(t,  x)  +  72  M4(t,  x), 

where  it  is  noted  that  since  A  depends  on  the  applied  magnetic  field,  it  is  not  homogeneous  throughout 
the  rod.  Hence,  the  magnetostriction  varies  along  x. 
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The  longitudinal  rod  displacements  u(t,  x )  are  computed  by  solving  equation  (22) 


d^u 

pAYt dx 


f 


du  , 


E  A  7t-(T  x)  +  cD  A 


d 2 


U 


( t .  a;')  —  E  A  Aft. 


d2u 

kL  u(t ,  L)  +  cL  —  (t,  X)  +  m/  -^-(f,  L) 


m- 


To  approximate  the  solution  to  this  equation,  a  Galerkin  discretization  in  x  is  used  to  reduce  the  system 
to  a  temporal  system  which  is  then  solved  with  finite  difference  approximations.  Details  regarding  the 
solution  method  used  are  provided  in  [11,  12],  Once  the  displacements  have  been  characterized,  the  strain 
is  computed  directly  using  the  equality 

/  \  9u 

£(M)  = 

and  the  corresponding  stresses  acting  on  the  rod  are  computed  directly  from  the  strain  using  equation  (21). 

The  complete  elastic-magnetic  coupled  model  is  summarized  in  Table  1.  It  should  be  noted  that  for 
implementation  purposes  it  is  convenient  to  assume  e  pa  A  for  computing  a  and  o7,  since  the  magnetostriction 
A  and  its  derivatives  with  respect  to  M  can  be  formulated  explicitly. 


6  Model  Validation 

In  this  section,  the  model  summarized  in  Table  1  is  applied  to  two  distinct  example  cases:  (i)  characteriza¬ 
tion  of  rod  magnetization  and  strain  output  by  a  magnetostrictive  transducer,  and  (ii)  characterization  of 
force  generated  by  a  magnetostrictive  rod  under  mechanically  blocked  (or,  equivalently,  electrically  short 
circuited)  conditions. 

The  first  example  demonstrates  the  model  accuracy  for  quantifying  the  quasistatic  magnetization  and 
strain  generated  by  a  magnetostrictive  material  when  used  in  a  prototypical  transducer.  The  performance 
of  the  model  is  evaluated  for  four  combinations  consisting  of  two  prestress  levels  and  two  drive  current 
levels.  It  is  shown  that  while  minor  changes  in  the  coupling  coefficient  a  and  magnetostriction  parameters 
7 i  are  necessary  to  fine  tune  fits  across  prestress  levels,  a  given  set  of  parameters  provides  accurate  fits 
across  drive  levels.  The  variation  in  a  across  prestress  levels  is  explained  by  the  fact  that  the  anhysteretic 
magnetization  is  stress  dependent,  and  that  the  shape  of  the  anhysteretic  has  a  strong  effect  on  the  shape 
of  the  prevailing  magnetization.  Regarding  the  magnetostriction  parameters  7 4-,  it  is  noted  that  the  domain 
configuration  does  vary  when  the  rod  is  unloaded  and  preloaded  again  for  operation  at  a  new  prestress 
value.  It  is  then  expected  that  these  domain  configuration  changes  be  accounted  for  by  variations  of 
parameters  7 4-. 

A  variety  of  transducer  applications  require  high  forces  of  the  kind  typically  associated  with  mechanically 
blocked  or  near-blocked  conditions.  Examples  of  such  applications  include  shaft  damping  in  inchworm- 
type  linear  motors  [35],  mitigation  of  seismic  vibrations  in  buildings,  and  deep-submersion  underwater 
communications.  This  motivates  the  evaluation  of  the  stress  calculations  provided  by  the  model.  Hence, 
in  the  second  example,  model  simulations  of  transducer  force  output  are  compared  with  experimental 
measurements  performed  under  mechanically  blocked  conditions  (s(t,  L)  =  0).  Experimental  tests  were 
performed  for  four  different  drive  current  levels,  at  a  constant  prestress.  It  is  shown  that  a  fixed  set  of 
parameters  provides  accurate  characterization  of  the  force  produced  by  a  magnetostrictive  rod. 

6.1  Example  1 

The  model  is  applied  to  a  Terfenol-D  transducer  with  the  configuration  illustrated  in  Figure  1  and  used  to 
characterize  the  quasistatic  magnetization  and  displacement  produced  by  the  transducer  in  response  to  an 
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Table  1:  Coupled  structural-magnetic  model  quantifying  the  magnetization,  strain,  and  stress. 


M(t,x)-(l-c)  Mirr{t ,  x)  +  c  [Man(t,  x)  +  Man(t,  x ) 

Structural  Model: 

A [M(t,  x)]  -  71  M2(t ,  x)  +  72  M4(t,  x) 

fL  02u  .  f1,  [  9m  .  d2u  1  dcj). 

Jo  PA-Qp(t-x)<l)(x)dx  =  ~  Jo  EA—{t,x)  +  cDA-Q-^{t,x)-EA\{t,x)  g^{x)dx 

,  /  n  9m  d2u 

-  kL  u[t,  L)  +  cL  —  (t,  L)  +  mL  L)  4>{L) 

.  .  9  m  . 

e(M)  = 

,  .  ,  ■  <)u ,  .  92m  ,  .  , .  , , 

cr(t,a;)  =  E  —(t,x)  +  cD  ^  A(t,  x)  +  <r0 


applied  current  I(t).  The  measured  output  from  the  transducer  during  operation  included  driving  voltage 
and  current,  voltage  induced  in  the  pickup  coil,  and  rod  displacement.  The  prestress  levels  were  —6.9  and 
—  10.35  MPa  (  —  1.0  and  —1.5  ksi),  and  the  drive  current  levels  were  4.0  and  16.0  A  zero-pk.  The  operating 
temperature  was  monitored  with  a  thermocouple  and  was  confined  to  the  range  21-28°C. 

The  driver  was  a  115  mm  (4.53  in)  long,  12.7  mm  (0.5  in)  diameter,  monolithic  Tbo.3Dyo.7Fe1. 92  rod 
manufactured  using  the  modified  Bridgman  process.  The  magnetic  field  excitation  H[i)  was  provided  by  a 
surrounding  solenoid,  while  an  innermost  one-layer  pickup  coil  was  used  to  collect  magnetic  induction  mea¬ 
surements.  The  1  Hz  excitation  frequency  was  provided  by  a  Tektronix  2642A  Personal  Fourier  Analyzer 
connected  to  a  Techron  7780  amplifier  operating  in  current  control  mode.  The  magnetic  induction  B{t) 
was  calculated  by  integration  of  the  pickup  coil  signal  Vpu.  Following  the  Faraday-Lenz  law  of  magnetic 
induction,  B(t)  =  —  1  /(Npu  Apu )  Vpu(r)  dr,  where  Apu  is  the  mean  cross  sectional  area  of  the  pickup  coil 
and  Npu  is  the  number  of  turns  in  the  pickup  coil.  The  magnetization  used  for  comparison  with  model 
simulations  was  computed  from  the  magnetic  constitutive  equation  M  =  (B/fj, 0)  —  H. 

Flux  closure  was  provided  by  1018  steel  end  caps  and  a  cylindrical  Alnico  V  permanent  magnet,  which 
was  demagnetized  to  obtain  unbiased  operation.  A  slider  and  a  Belleville  compression  washer  completed 
the  magnetic  circuit. 

The  compressive  preloading  necessary  to  avoid  tensile  stresses  on  the  moderately  brittle  Terfenol-D  rod 
was  provided  by  a  steel  bolt  located  within  the  transducer  base.  The  bolt  pushed  the  rod  against  the 
compression  washer,  whose  stiffness  kmps  was  calculated  from  the  linear  region  of  its  force-displacement 
characteristic  curve.  A  PCB  208A13  force  transducer  located  between  the  prestress  bolt  and  the  Terfenol-D 
rod  was  used  during  preloading  to  obtain  the  desired  preload  value. 

The  external  load  was  a  mass  threaded  onto  the  displacement  plunger  which  weighed  m ^  =  0.5  Kg. 
The  load  displacement  (i.e.,  the  transducer  output)  was  measured  with  a  Lucas  Schaevitz  LVM-110  linear 
variable  differential  transducer  equipped  with  a  025-MHR  probe.  The  corresponding  bulk  strain,  used  for 
comparison  with  model  simulations,  was  computed  by  dividing  this  displacement  by  L. 

For  identification  of  the  $  filter,  a  series  of  flux  measurements  was  conducted  with  a  Hall  effect  probe  (F. 
W.  Bell  9500  series),  located  within  the  transducer  to  allow  characterization  of  end  effects,  demagnetizing 
factors,  and  magnetic  circuit  nonlinearities.  A  DC  current  was  used  to  excite  the  solenoid,  and  the  flux 
density  Bsur(x )  was  measured  adjacent  to  the  surface  of  the  rod  at  locations  5  mm  apart  over  a  longitudinal 
line.  Since  in  air  B  =  H  in  the  CGS  units  system,  the  flux  density  is  equal  to  the  magnetic  field, 
BSur{x)  =  Hsur(x).  The  correction  functional  was  then  computed  from  T(a:)  =  Hsur{x )  /I.  Additional 
details  regarding  the  experimental  analysis  are  provided  in  [11,  12], 

The  magnetization  and  strain  calculations  measured  from  the  experimental  transducer  are  shown  in 
Figure  5.  It  should  be  noted  that  the  reversal  of  strain  magnitude  observed  across  drive  levels  when  the 
prestress  is  changed  from  —6.9  MPa  to  —10.35  MPa  is  consistent  with  the  data  reported  in  [36]. 

The  performance  of  the  model  is  illustrated  in  Figures  6  and  7.  The  model  provides  in  all  cases  very 
accurate  fits  of  both  the  measured  rod  magnetization  and  transducer  strain  output,  ft  should  be  noted 
that  the  model  does  not  provide  a  mechanism  to  account  for  the  contraction  observed  in  the  magnetization 
data  at  low  field  levels,  which  explains  the  discrepancy  observed  in  the  4.0  A  cases.  It  is  expected  that  a 
higher  order  truncation  of  the  magnetostriction  expansion  should  translate  into  improved  strain  fits. 

6.2  Example  2 

The  model  is  now  used  to  characterize  the  force  output  by  a  19.05  mm  (0.75  in)  long,  8.89  mm  (0.35  in) 
diameter,  monolithic  Tbo.3Dyo.7Fe1 .92  FSZM  rod.  A  testing  clamp  was  devised  to  maintain  mechanically 
blocked  conditions  during  operation,  which  consisted  of  a  stainless  steel  base  where  the  rod  and  magnetic 
circuit  assembly  were  located  and  a  stainless  steel  cap  which  housed  a  PCB  A13  force  transducer.  The  cap 
was  attached  to  the  base  via  stainless  steel  bolts,  which  were  tightened  until  the  desired  prestress  value 
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1 


(a) 


(b) 

Figure  5:  Experimental  data  collected  from  a  Terfenol-D  transducer,  at  the  prestress  values 
—  6.9/  — 10.35  MPa  and  the  drive  current  levels  4.0/16.0  A  zero-pk.  (a)  Magnetization  and  (b)  output 
strain. 
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Input  Current  4.0  A  zero-pk 
Prestress  —6.9  MPa 


Input  Current  4.0  A  zero-pk 
Prestress  —10.35  MPa 


H  [kA/m] 


H  [kA/m] 


Input  Current  16.0  A  zero-pk 
Prestress  —6.9  MPa 


Input  Current  16.0  A  zero-pk 
Prestress  —10.35  MPa 


H  [kA/m]  H  [kA/m] 


Figure  6:  Magnetization  results:  comparison  of  model  simulations  ( - )  with  experimental 

data  (—•—).  Constant  values  of  the  model  parameters  were  a  —  7000  A/m,  k  —  7000  A/m,  c  = 
0.2,  a  =  0.045,  E  =  30  GPa,  p  =  9250  Kg/m3,  7l  =  2.95xl0“15  m2/A2,  y2  =  -6xl0“28  m4/A4, 
£  =  2.45  X  104  Pa,  cd  —  lx  106  Ns/m,  cj=lx  103  Ns/m,  =  2.5  X  106  N/m,  m./,  —  0.5  Kg, 
L  =  0.115  m,  D  =  0.0127  m. 
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Input  Current  4.0  A  zero-pk 
Prestress  —6.9  MPa 


Input  Current  4.0  A  zero-pk 
Prestress  —10.35  MPa 


Input  Current  16.0  A  zero-pk 
Prestress  —6.9  MPa 


Input  Current  16.0  A  zero-pk 
Prestress  —10.35  MPa 


Figure  7:  Strain  results:  comparison  of  model  simulations  ( - )  with  experimental  data  ( - ). 

Constant  values  of  the  model  parameters  were  same  as  in  Figure  6.6,  except  a  =  0.065,  71  = 
3.12  X  10-15  m2/A2. 
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was  read  from  the  force  transducer. 

The  magnetic  field  H{t)  was  supplied  by  a  400-turn  solenoid,  which  was  calibrated  using  a  Hall  effect 
probe  in  the  fashion  described  in  Example  f.  The  magnetic  circuit  was  completed  by  two  magnetic  steel 
end  pieces  and  a  surrounding  magnetic  steel  cylinder. 

The  tests  were  performed  at  a  prestress  value  of  —6.9  MPa  for  four  drive  current  levels:  0.5,  f.0,  1.5, 
and  2.0  A  zero-pk.  The  excitation  signal  was  1  Hz  in  all  cases.  The  operating  temperature  remained  in  all 
cases  between  2f-30°C.  The  force  versus  magnetic  field  measurements  are  shown  in  Figure  8. 


Figure  8:  Force  versus  magnetic  field  generated  by  the  Terfenol-D  rod  under  blocked  conditions,  for  the 
drive  levels  0.5,  f.0,  1.5,  and  2.0  A  zero-pk.  The  prestress  was  <7o  =  —6.9  MPa. 

Although  in  theory  the  clamp’s  stiffness  should  be  infinity,  the  manufacturer’s  specifications  for  the 
force  transducer  were  commensurate  with  that  of  the  Terfenol-D  rod  itself.  To  model  this  condition,  the 
specified  stiffness  value  was  assigned  to  the  boundary  spring,  kmps  =  f.75  X  f08  N/m.  In  a  similar  fashion, 
the  dynamic  mass  was  weighed  and  assigned  to  the  corresponding  model  mass  load,  rn /.  =  5  Kg. 

ft  should  be  noted  that  the  elastic  modulus  under  blocked  conditions  (or  equivalently,  under  constant 
magnetic  induction  conditions)  represents  the  intrinsic  material  stiffness  and  is  the  largest  value  that  the 
material  can  achieve  for  a  given  operating  condition.  The  value  of  the  elastic  modulus  used  in  these 
simulations  was  E  =  60  GPa,  in  accordance  with  the  experimental  measurements  reported  in  [34]. 

Model  simulations  and  experimental  measurements  are  compared  in  Figure  9.  The  model  performance 
is  extremely  satisfactory  at  all  drive  levels,  both  in  quantifying  the  amplitude  of  the  force  and  the  hysteresis 
present  in  the  data,  ft  is  emphasized  that  the  same  set  of  values  for  model  parameters  was  used  in  all 
cases. 
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Input  Current  0.5  A  Input  Current  1.0  A 


Input  Current  1.5  A  Input  Current  2.0  A 


Figure  9:  Comparison  of  model  simulations  ( - )  with  experimental  data  (—  ■  — )  at  four  drive  levels 

and  prestress  of  oo  =  —6.9  MPa.  The  constant  values  of  the  model  parameters  were  a  =  4000  A/m, 
k  =  5000  A/m,  c  =  0.4,  a  =  0.001,  E  =  60  GPa,  p  =  9250  Kg/m3,  7l  =  7  X  10“16  m2/A2,  72  = 
—  1  X  10-29  m4/A4,  £  =  2.45  X  104  Pa,  cD  =  1  X  105  Ns/m,  cL  =  1  X  103  Ns/m,  —  1.75  X  10s  N/m, 
m.L  =  5  Kg,  L  =  0.019  m,  D  =  0.009  m. 
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7  Concluding  Remarks 


A  coupled  magnetomechanical  model  for  the  magnetization  and  strain  behavior  of  magnetostrictive  trans¬ 
ducers  in  response  to  applied  currents  has  been  presented.  The  model  includes  the  nonlinearities  and  hys¬ 
teresis  present  in  the  magnetic  response  of  magnetostrictive  materials,  and  the  linear  elastic  effects  which 
are  characteristic  of  high-signal  magnetostrictive  transduction.  The  magnetic  and  structural  regimes  are 
coupled  through  the  magnetization,  which  varies  in  response  to  both  the  externally  applied  magnetic  field 
and  the  stress  field  which  arises  as  the  transducer  actuates  upon  an  external  load. 

The  model  was  constructed  in  three  steps.  In  the  first,  the  mean  field  theory  of  Jiles  and  Atherton 
was  used  as  a  basis  for  quantifying  the  relation  between  the  current  input  to  the  transducer  and  the 
magnetization.  The  magnetization  arising  from  variations  in  stress  was  also  addressed,  by  considering 
a  ‘law  of  approach’  to  the  anhysteretic  magnetization.  The  resulting  magnetization  model  provides  a 
representation  of  the  hysteresis  and  saturation  effects  taking  place  when  domain  walls  attach  to  and 
unattach  from  pinning  sites  in  the  material.  Pinning  sites  provide  a  mechanism  for  restraining  the  motion 
of  domain  walls  as  the  field  is  cycled,  and  hence  for  magnetization  hysteresis.  The  stress  acting  on  the 
material  produces  domain  wall  unpinning  and  hence  causes  the  magnetization  to  approach  the  anhysteretic 
state.  In  the  second  step,  the  magnetostriction  due  to  the  rotation  of  magnetic  moments  was  quantified 
by  means  of  a  quartic  model  formulated  in  terms  of  the  magnetization.  Finally,  force  balancing  was  posed 
in  the  form  of  a  PDE  equation  which  includes  the  magnetostriction,  system  compliance,  internal  damping, 
and  boundary  conditions  given  by  the  transducer  design.  The  solution  to  this  PDE  provides  the  material 
displacements  and  forces  produced  by  the  transducer. 

The  examples  demonstrated  the  accuracy  of  the  model  in  two  cases  typical  of  magnetostrictive  transducer 
applications.  In  the  first,  the  magnetization  and  strain  output  by  a  typical  actuator  were  characterized 
at  two  drive  levels  and  two  prestress  levels.  A  near  constant  set  of  parameters  was  used  across  drive 
levels  and  prestresses.  In  the  second  example,  the  force  output  by  a  magnetostrictive  rod  operated  under 
mechanically  blocked  conditions  was  characterized  for  four  drive  levels  and  fixed  prestress.  A  constant  set 
of  model  parameters  provided  accurate  simulations  across  drive  levels. 
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